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Two-Loop Sudakov Form Factor in a Theory with Mass Gap 



in 
o 
o 

(N 
< 



(N 
> 

00 ■ 

o ■ 
^t- : 
o . 

o ; 
^ : 
i 

^ ■ 



X 



Bernd Feucht, 1,2 Johann H. Kuhn, 2 ' 3 Alexander A. Penin, 2 ' 4 and Vladimir A. Smirnov 5,6 

1 Bernd Jantzen in later publications 
2 Institut fur Theoretische Teilchenphysik, Universitat Karlsruhe, 76128 Karlsruhe, Germany 
Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106 USA 
** Institute for Nuclear Research, Russian Academy of Sciences, 117312 Moscow, Russia 

5 Institute for Nuclear Physics, Moscow State University, 119992 Moscow, Russia 
6 II. Institut fur Theoretische Physik, Universitat Hamburg, 22761 Hamburg, Germany 

The two-loop Sudakov form factor is computed in a 17(1) model with a massive gauge boson and 
a U(l) x U(l) model with mass gap. We analyze the result in the context of hard and infrared 
evolution equations and establish a matching procedure which relates the theories with and without 
mass gap setting the stage for the complete calculation of the dominant two-loop corrections to 
electroweak processes at high energy. 
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Since the pioneering works by Sudakov and Jackiw 
the high energy asymptotics of the electromagnetic 
form factor has been the subject of numerous investiga- 
tions. The problem is relevant for a wide class of phe- 
nomenological applications from Drell-Yan processes to 
deep inelastic scattering. Recently a new wave of inter- 
est to the Sudakov asymptotic regime has been risen in 
connection with higher-order corrections to electroweak 
processes at high energies 0, H 0, 0, 0, S S E3, - Ex- 
perimental and theoretical studies of electroweak inter- 
actions have traditionally explored the range from very 
low energies, e.g. through parity violation in atoms, up 
to energies comparable to the masses of the W- and Z- 
bosons, e.g. at the LEP or the Tevatron. The advent 
of multi-TeV colliders like the LHC during the present 
decade or a future linear electron-positron collider will 
give access to a completely new energy domain. Once the 
characteristic energies y/s are far larger than the masses 
of the W- and Z-bosons, Mw,z, exclusive reactions like 
electron-positron (or quark-antiquark) annihilation into 
a pair of fermions or gauge bosons will receive virtual 
corrections enhanced by powers of the large electroweak 
logarithm ln(s/M^ z ). The leading double-logarithmic 
corrections may well amount to ten or even twenty per- 
cent in one-loop approximation and reach a few percent 
in two- loop approximation. Moreover, in the TeV re- 
gion, the subleading logarithms turn out to be equally 
important [(J [^j and a percent accuracy of the theo- 
retical estimates for the cross sections necessary for the 
search of new physics beyond the standard model can be 
guaranteed only by including all the logarithmic two- loop 
corrections. 

The calculation of the two-loop electroweak corrections 
even in the high energy limit is an extremely challenging 
theoretical problem. It is complicated in particular by 
the presence of the mass gap and mixing in the gauge 
sector. However, the logarithmic corrections are quite 



insensitive to fine details of the spontaneous symmetry 
breaking. The calculation of the leading logarithmic (LL) 
electroweak corrections can be performed using the fields 
of the unbroken symmetry phase and how the infrared 
singular virtual photon contribution can be separated 
within the infrared evolution equation approach |5( . This 
scheme has been extended to the next-to-leading (NLL) 
and next-to-next-to-leading logarithmic (N 2 LL) approx- 
imation in Refs. [rjHo|. 

In the study of Sudakov corrections the analysis of the 
form factor plays a special role since it is the simplest 
quantity which includes the complete information about 
the universal collinear logarithms 12] directly applicable 
to a process with an arbitrary number of fermions. In 
this Letter we formulate a general matching procedure 
which relates the logarithmic corrections in the theories 
with and without mass gap by combining the hard and 
infrared evolution equation approach with the explicit 
two-loop results for the form factor in an Abelian gauge 
model. This reduces the calculation of the dominant two- 
loop corrections to electroweak processes at high energy 
to a single-mass problem without mixing. 

The structure of the Letter is as follows. First, we 
present the explicit two-loop results for the form factor 
in a U(l) model with a massive gauge boson. Then we 
introduce the evolution equations, compute the two-loop 
corrections to the form factor in a U(l) x U(l) model 
with mass gap, and establish the matching procedure. 
Finally we outline how the approach can be applied to 
the calculation of the two- loop electroweak corrections to 
neutral current four-fermion processes. 

The vector form factor T determines the fermion scat- 
tering amplitude in an external Abelian field. It is 
a function of the Euclidean momentum transfer Q 2 — 
~(pi ~ P2) 2 where pi t 2 is the incoming/outgoing fermion 
momentum and we consider on-shell massless fermions, 
p\ = p\= 0. Let us write the perturbative expansion for 



2 



0.010 



0.005 



0.000 




500 



1000 



1500 2000 

Q [GeV] 



2500 



3000 



0.010 



0.005 



0.000 




500 



1000 



1500 2000 

Q [GeV] 



2500 



3000 



FIG. 1: The two- loop correction to the form factor J- a (M, Q) 
in LL (including a 2 £ 4 ), NLL (including a 2 £ 3 ), N 2 LL (includ- 
ing a 2 C 2 ), N 3 LL (including a 2 C 1 ) approximations and the 
complete two-loop correction as functions of the momentum 
transfer for M = 80 GeV, a /{Ant) = 3 ■ 10 -3 . 



the form factor as F a {M, Q) = (^)" ] (n) T B where 
Tb corresponds to the Born approximation, /'°) = 1. 
For the U{1) model with a gauge boson of mass M in the 
Sudakov limit M/Q — ► the one-loop correction is well 
known 
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where £ = In (Q 2 /M 2 ) and all the power-suppressed 
terms are neglected. For the two-loop term we find 
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where C(3) = 1.202057 ... and Li 4 (§) = 0.517479 . . . are 
the values of the Riemann's ^-function and the polylog- 
arithm, respectively. In Eq. (0 we do not include the 
contribution due to the virtual fermion loop computed 
in 0| ■ For the calculation of the leading power behavior 
of the two-loop on-shell vertex diagrams with two massive 
propagators in the Sudakov limit we used the expansion 
by regions approach (for the application to the Su- 
dakov form factor see also Q). The method is based on 
the factorization of the contributions of the dynamical 
modes characteristic for the Sudakov limit 0] in dimen- 
sional regularization. Our result for the contribution of 
the hard modes agrees with the dimensionally regularized 
massless result of Ref. . Details of our calculation will 
be published elsewhere. 

In Fig. n the numerical results for the two-loop cor- 
rection to the form factor in the different logarithmic 



FIG. 2: The same as Fig. Rafter changing the argument of 
the logarithm. 



approximations are plotted as functions of the momen- 
tum for the values of M and a typical for electroweak 
interactions. The two-loop logarithmic terms have a 
sign-alternating structure resulting in significant cancel- 
lations. In the region of a few TeV the form factor does 
not reach the double-logarithmic asymptotics. The quar- 
tic, cubic and quadratic logarithms are comparable in 
magnitude and dominate the two-loop corrections. Then 
the logarithmic expansion starts to converge and, after 
including the linear-logarithmic contribution, provides a 
very accurate approximation of the total two-loop cor- 
rection. Such a behavior is typical for the Sudakov limit 
and holds for the non-Abelian corrections as well j 1 OL 1 1 1| . 
Note that by rescaling M — > e 3 / 4 M in the argument of 
the logarithm the NLL contribution can be made to dis- 
appear. That improves significantly the convergence of 
the logarithmic expansion and prevents the strong can- 
cellation between the logarithmic terms (see Fig.|2J). Still, 
the N 3 LL contribution is a must for the quantitative ap- 
proximation. 

The asymptotic dependence of the form factor on Q is 
governed by the linear hard evolution equation [lri | . As 
a consequence, the logarithmic corrections exponentiate. 
For the purely Abelian contribution the exponent has a 
particularly simple form 
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The double-logarithmic term in the exponent is pro- 
tected against the Abelian multiloop corrections by the 
properties of the light-cone Wilson loop 17] . Our two- 
loop result determines the next-to-next-to-next-to-next- 
to-leading logarithmic (N 4 LL) approximation of the form 
factor which includes the a n C m corrections with m = 
2n — 4, . . . , In to all orders in a. 
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Let us now turn to the second example with two 
Abelian gauge bosons of the masses A and M, A <C M, 
and couplings a' and a, respectively. We can introduce 
the infrared evolution equation which governs the depen- 
dence of the form factor T(X, M, Q) on A The virtual 
corrections become divergent in the limit A — > 0. Accord- 
ing to the Kinoshita-Lee-Nauenberg theorem these 
divergences are cancelled against the ones of the correc- 
tions due to the emission of real light gauge bosons of 
vanishing energy and/or collinear to one of the on-shell 
fermion lines. The singular behavior of the form factor 
must be the same in the full U a >(l) x U a (l) theory and 
the effective U a t(l) model with only the light gauge bo- 
son. Thus for A <C M <§C Q the solution of the infrared 
evolution equation is given by the exponent of Eq. © 
with M, a replaced by A, a', and the form factor can be 
written in a factorized form 

F{\ M, Q) = F(M, Q)F a , (A, Q) + 0(X/M) , (4) 

where the function F{M, Q) depends on a and a' and 
incorporates all the logarithms of the form In (Q 2 /M 2 ). 
It can be obtained directly by calculating the ratio 
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Since the function F(M, Q) does not depend on the in- 
frared regularization, we compute the ratio in Eq. (J5J 
with A = using dimensional regularization for the in- 
frared divergences. The method of calculation of the 
two-loop diagrams with both massive and massless gauge 
bosons is similar to the purely massive case. We obtain 
the two-parameter perturbative expansion F(M, Q) = 
E„, m ^£r/ (n ' m) > where /C°.°> = 1, /K°> = 0, 
y(o,m) _ an d the two-loop interference term reads 
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The numerical structure of the corrections to F(M, Q) is 
very similar to the one of F a {M, Q) (see Fig.[Q. 

In the equal mass case, A = M, we have an additional 
reparameterization symmetry, and the form factor is de- 
termined by Eq. J2J) with the effective coupling a — a' +a 
so that F(M,M,Q) = F & (M, Q). We can now write 
down the matching relation 



T{M, M, Q) = C(M, Q)F(M, Q)T a , (M, Q) 



(7) 



where the matching coefficient C (M, Q) represents the ef- 
fect of the power-suppressed terms neglected in Eq. 
By combining the explicit results for J- a i(M,Q) and 
F(M, Q) we find the two-loop matching coefficient 
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Eq. ||SJ| does not contain logarithmic terms, and up to the 
N 3 LL accuracy the product F(M,Q)J 7 a /(X,Q) continu- 
ously approaches T{M, M, Q) as A goes to M. There- 
fore, to get all the logarithms of the heavy gauge boson 
mass in two-loop approximation for the theory with mass 
gap, it is sufficient to divide the form factor T a {M,Q) 
of the symmetric phase by the form factor T a > (A, Q) of 
the effective U a i(l) theory taken at the symmetric point 
A = M. Thus we have reduced the calculation in the 
theory with mass gap to the one in the symmetric theory 
with a single mass parameter. Note that the absence of 
the linear-logarithmic term in Eq. JSJ is an exceptional 
feature of the Abelian corrections. The general analy- 
sis of the evolution equation 0] shows the terms ne- 
glected in Eq. (@J to contribute starting from the N 3 LL 
approximation. This implies the absence of the second 
and higher powers of the logarithm in the matching co- 
efficient of Eq. Q , irrespectively of the gauge group and 
the mass generation mechanism. Moreover, in the ap- 
proximately equal mass case, (M — X)/M = 5 <C 1, one 
can compute the form factor as an expansion around the 
equal mass result. Up to N 2 LL accuracy only the leading 
term of Eq. contributes and the expansion takes the 
form 



xT & (M,Q) + 0{Sa'a£) . (9) 



Let us show how the above procedure applies to the calcu- 
lation of two-loop electroweak corrections. To be specific, 
we consider a four- fermion neutral current process, which 
is of primary phenomenological importance, with light 
fermions. The four- fermion amplitude can be decom- 
posed into (the square of) the form factor and a reduced 
amplitude ja, 0] • The latter carries all the Lorentz and 
isospin indices and does not contain collinear logarithms 
in perturbative expansion. The logarithmic corrections 
to the reduced amplitude are obtained by solving a renor- 
malization group like equation ^^|. The corresponding 
two-loop anomalous dimensions can be extracted from 
the existing massless QCD calculations HI (see USED). 
Thus, the problem of the calculation of the two-loop elec- 
troweak logarithms in the four-fermion processes reduces 
to the analysis of the form factor. 

In Ref. [l(J by analyzing the hard evolution equation 
it has been found that the two-loop electroweak correc- 
tions up to the next-to-next-to-leading (quadratic) log- 
arithms are not sensitive to the structure of the theory 
at the electroweak symmetry breaking scale. The predic- 
tion of Ref. 0, ^3 for the two-loop logarithmic correc- 
tions fully agrees with the available explicit results for 
the light fermion contribution [xil | and the Abelian con- 
tribution obtained in this Letter. The only trace of the 
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Higgs mechanism of the gauge boson mass generation in 
N 2 LL approximation is the Z — W boson mass splitting 
which can be systematically taken into account within an 
expansion around the equal mass approximation similar 
to Eq. @. Thus, the calculation of the two-loop elec- 
troweak corrections up to the quadratic logarithms can 
be performed in two steps outlined above: (i) the cor- 
rections are evaluated using the fields of unbroken sym- 
metry phase with all the gauge bosons of the same mass 
M s» Mz,w introduced by hand; (ii) the QED contribu- 
tion with an auxiliary photon mass M is factorized as in 
Eq. J7J leaving the pure electroweak logarithms. The sep- 
arated virtual QED corrections accompanied by the real 
photon radiation in the limit of vanishing photon mass 
result in the universal infrared safe factor independent of 
M z ,w. 

By contrast, the N 3 LL approximation is sensitive to 
fine details of the gauge boson mass generation and the 
coefficient of the linear two-loop electroweak logarithm 
depends e.g. on the Higgs boson mass. For the full 
calculation of this coefficient one has to use the true 
mass eigenstates of the standard model. Our result, 
Eqs. EJ), is an example of such a calculation when 
applied to the two-loop diagrams with photon and Z bo- 
son exchanges. We can, however, make a reasonable ap- 
proximation which dramatically simplifies the analysis. 
Namely, consider a simplified model with a Higgs bo- 
son of zero hypercharge. Then the mixing is absent and 
the hypercharge gauge boson remains massless. The in- 
terference diagrams including the heavy SUl{2) and the 
light hypercharge U(l) gauge bosons are identical with 
the ones of the purely Abelian model discussed in this 
Letter, where the above two-step procedure can be ap- 
plied to get all the two-loop logarithms including the lin- 
ear term. In the standard model the mixing of the gauge 
bosons results in a linear-logarithmic contribution, which 
is not accounted for within this procedure. It is, however, 
suppressed by a small factor sin 2 9w ~ 0.2, with 8w be- 
ing the Weinberg angle. Therefore, the above simplified 
model gives an estimate of the coefficient in front of the 
linear electroweak logarithm with 20% accuracy. From 
the numerical result of Fig. ^ which represents the typ- 
ical structure of the two-loop corrections, we see that a 
20% error in this coefficient leads to an uncertainty com- 
parable to the nonlogarithmic contribution and is prac- 
tically negligible. Thus we are able to get an accurate 
estimate of the two-loop correction, which is sufficient 
for practical applications to the future collider physics, 
by performing the calculations in the model without mass 
gap and mixing of the gauge bosons. The last ingredient 
necessary to complete the calculation of the dominant 
two-loop electroweak corrections is the generalization of 
Eq. (J2J) up to the linear-logarithmic term to the pure 
SUl (2) gauge model with the Higgs mechanism of mass 
generation. Note that up to the quadratic logarithm the 
two-loop corrections are predicted by the evolution equa- 



tion 0. 

To conclude, we have obtained the complete results 
for the two-loop corrections to the vector form factor in 
the Sudakov limit in Abelian theories with one massive 
gauge boson or with two gauge bosons of essentially dif- 
ferent masses. The results are in full agreement with 
the predictions of the evolution equation approach. We 
have formulated a systematic procedure of factorizing 
the infrared singular virtual corrections and reducing the 
calculation in the theory with mass gap to the single- 
mass problem. The analysis can directly be generalized 
to the standard model with spontaneous breaking of the 
SU L {2) xU(l) theory to low-energy QED. This solves 
the principal problems of the calculation of the domi- 
nant two-loop electroweak corrections to the neutral cur- 
rent four-fermion processes which are mandatory for the 
high-precision physics at the LHC and the next genera- 
tion of linear colliders. 
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